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Abstract. Let 9Jl„.m be the space of real n x m matrices which 
can be identified with the Euchdean space R"™. We introduce con- 
tinuous wavelet transforms on ^n,m with a multivalued scaling pa- 
rameter represented by a positive definite symmetric matrix. These 
transforms agree with the polar decomposition on SJln.m and coin- 
cide with classical ones in the rank-one case m — 1. We prove an 
analog of Calderon's reproducing formula for i^-functions and ob- 
tain explicit inversion formulas for the Riesz potentials and Radon 
transforms on dJln,m- We also introduce continuous ridgelet trans- 
forms associated to matrix planes in ^IIln^m■ An inversion formula 
for these transforms follows from that for the Radon transform. 



It is known that diverse wavelet-like transforms can be generated by 
operators of fractional integration and used to invert these operators. 
On the other hand, numerous problems in integral geometry, for in- 
stance, reconstruction of functions from their integrals over planes in 
M", reduce to inversion of fractional integrals. The following example 
illustrates these statements and explains how wavelet transforms arise 
in the context of integral-geometrical problems; see also )Ru2j for the 
more detailed exposition. 
Consider the Riesz potential 
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7]-"/22"r(Q;/2) 
r{{n-a)/2) 
and replace the kernel \x — yj^^" by the integral 

\x — y\\ da 



oo 



;i.2) \x - yr^^ = / w 



where w(-) is good enough and c^,^ = w{s) " ^ds ^ 0. Chang- 
ing the order of integration, we obtain 

-1 °° 

(1-3) a /)(^) = y ^T^^«' 



where 

(1.4) (W,/)(x) = f{y)w (^^) 

If w obeys some cancellation conditions, then ()1.4|1 represents the clas- 
sical continuous wavelet transform with the scaling parameter a > 
|Daj . |F JWj ■ |Hoj . If we start with a fractional integral different from 
/"/, say, with the Bessel potential or whatever (see |Rulj . Section 
10.7), we arrive at a wavelet transform, which differs from ()1.4|) . 

Since the Fourier transform of /"/ is \y\~°'{J^ f){y) in a certain sense, 
then, formally, (/")^^ = and it is natural to expect that the inverse 
operator (/")~^ can be represented in the form ()1.H|1 with a replaced 
by —a, namely. 



oo 



(1.5) (n-V = / ^c^a. 



da,w being a normalizing factor. In particular, for a = 0, 

oo 

(1.6) / = c/o,^ f^^^da. 



The equality (|1.6p is a modification of Calderon's reproducing formula 
|Calj ■ |FJWj, |iRu3j . |Ru7j . If w is normalized and has the form w = 
u*v, then (jl.fjj] turns into the classical Calderon identity 

oo 

(1.7) f=[^^^^da 



where Ua{x) = a "'u{x/a), Va{x) = a "'v{x/a). 
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Of course, this argument is purely heuristic and formulas ()1.5|l - ()1.7j) 
require justification in the framework of a suitable class of functions / 
under certain cancellation conditions for the wavelet function w. 

What is the connection between this argument and the Radon trans- 
form in integral geometry? Suppose that T is the manifold of all k- 
dimensional planes r in M". For functions / : M" C and (/? : T — C, 
the Radon transform and its dual are defined by 



(1-8) /(r) = j fix), = J ^(r), 

respectively, and obey the Fuglede equality 

(1.9) {fy = cl'f, c = const, 

see |Fuj . |Helj . |Ru4j for details. Combining ()1.9p with p.5|) . we obtain 
an inversion formula for / in the "wavelet form" 

oo 

(1.10) / = const j -^ri^da, 



where = limg^o in a certain sense. 

This formalism can be applied in a wider context, when, instead of 
the Euclidean distance |a; — ?/| between two points, one deals with the 
distance |a; — r| between the point a; G and the /c-dimensional plane 
r C M". Starting with the intertwining operator 

(1.11) (P"/)(r) = /"/(x)|a;-r|"+'=-"c/x 

7n-fc(a) J 

which is called the generalized Semyanistyi fractional integral (cf. jBi] 
for k = n — 1), one arrives at the corresponding wavelet-like trans- 
form, recently called the continuous k-plane ridgelet transform; see 
|Caj . |Ru5j ■ and references therein. These transforms have proved to 
be useful in applications |Caj . |Doj . |Murj . and are of independent the- 
oretical interest. 

A common feature of these examples is that the scaling parameter 
a > is one-dimensional no matter what the dimension of the ambient 
space is. The situation changes drastically if we replace n by nm, 
and regard M"™" as the space of n x m real matrices x = (xij). Then a 
similar procedure, starting with the properly defined Riesz potential, 
yields a new wavelet-like transform which is applicable to functions of 
matrix argument and depends on the matrix-valued scaling parameter. 
This parameter is represented by a positive definite symmetric matrix 
of size m X m. Apart from extra flexibility that might be useful in 
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applications, such transforms have a rich theory which rehes on diverse 
higher rank phenomena. 

In the present paper, we focus on the theory of the new wavelet 
transforms mentioned above, and give some applications. The paper 
is organized as follows. Section 2 contains necessary prerequisites. We 
fix our notation and recall basic facts related to Riesz potentials and 
Radon transforms on the space of rectangular matrices. In Section 
3, we introduce continuous wavelet transforms for functions of ma- 
trix argument and prove the corresponding reproducing formula of the 
Calderon type. In Section 4, we show how wavelet transforms can 
be used for inversion of Riesz potentials on matrix spaces. Unlike the 
rank-one case m = 1, for which numerous inversion formulas are known 
[Rul', [S KMj . the corresponding higher rank problem is very difficult; 
see |iQR2] . jRu6j for the discussion. Wavelet transforms prove to be a 
convenient tool to resolve this problem in the L^-case. In Section 5, 
we apply our wavelet transforms to inversion of the Radon transform 
associated to the so-called matrix /c-planes. These Radon transforms 
were studied in detail in jPe], |URlj . and |UR2j . where it was shown 
that the inversion problem for them has the same difficulties as for 
the Riesz potentials on the space of rectangular matrices. In Section 
5, we introduce continuous ridgelet transforms of functions of matrix 
argument, generalizing those in [Calj and 'Ru5j , and prove a reproduc- 
ing formula for these transforms. This result is a consequence of the 
inversion formula for the Radon transform. 

Acknowledgements. The work on the paper was started when B. 
Rubin was visiting Department of Mathematics at the Louisiana State 
University in April 2004. He is deeply grateful to his colleague. Prof. 
Gestur Olafsson for the hospitality. 

2. Preliminaries 

In this section, we fix our notation and recall some basic facts, that will 
be used throughout the paper. The main references are |Muj . [0R2j . 
p. 

2.1. Notation and some auxiliary facts. Let dJln,m be the space of 
real matrices x = (xij) having n rows and m columns. We identify 
Mn^m with the real Euclidean space M"™ and set dx = YYi=i YYjLi 
for the Lebesgue measure on Mn,m- In the following, x' denotes the 
transpose of x, Im is the identity m x m matrix, stands for zero 
entries. Given a square matrix a, we denote by tr(a) the trace of a, 
and by |a| the absolute value of the determinant of a, respectively. We 
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hope the reader will not confuse \a\ with the similar notation for the 
absolute value of a number because the meaning of a will be clear each 
time from the context. For x G 97ln,m, n > m, we set 

(2.1) \x\m = detix'xY^^ = \x'x\^^\ 

If m = 1, this is the usual Euclidean norm in M". If m > 1, then \x\m 
is the volume of the parallelepiped spanned by the column-vectors of 
the matrix x, cf. jHl P- 251]. 

Let Vm be the cone of positive definite symmetric matrices r = 
{rij)mxm with the elementary volume dr = Y[i<j ^''^ijy '^rn be 

the closure of Vm, that is the set of all positive semi-definite m x m 
matrices. We write r > if r G Vm, and r > if r G Vm, respectively. 
Given si and S2 in Vm, we write Si > $2 for si — S2 G Vm- If a G Vm 
and h G Vm, then f{s)ds denotes the integral over the compact set 

{a,b) = {s : s - a e Vm, b - s e Vm}- 

The group G = GL{m, M) of real non-singular m x m matrices g acts 
transitively on Vm by the rule r — > grg'. The corresponding G-invariant 
measure is 

(2.2) d^r = Irl-'^dr, |r| = det(r), d={m + l)/2 

m p- 18]. 

A function wq on Vm is called symmetric if 

(2.3) Wo{s'/\s^/^) = Wo{r'/hr^/^), Vr, s G Vm- 

For example, any function of the form Wo{r) = Wi(tr(r)) is sym- 
metric. Instead of the trace, one can take any function of the form 
ip{ai, . . . , am), where ai, . . . ,am are elementary symmetric functions of 
the eigenvalues Ai, . . . , Xm of r. This follows from the general fact that 
if A and B are nonsingular square matrices then AB and BA have the 
same eigenvalues; see, e.g., |Mul pp. 584, 585]. 

We use a standard notation 0{n) for the group of real orthogonal 
n X n matrices; SO{n) = {7 G 0{n) : det(7) = 1}. The corresponding 
invariant measures on 0{n) and SO{n) are normalized to be of total 
mass 1. The Lebesgue space = L^(9Jl„,m) and the Schwartz space 
S = §{Tln^m) are identified with respective spaces on M"'". We denote 
by CciMn^m) the space of compactly supported continuous functions on 

■-'■'^n,m- 

The Fourier transform of a function / G L^{Mn^m) is defined by 

(2.4) (^/)(l/)= / exp{tr{iy'x))f{x)dx, y E Tln,m - 
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This is the usual Fourier transform on R"™ so that the relevant Parseval 
formula reads 

(2.5) (^/,^¥,) = (27r)"'"(/,v^), 

where 



if,^)= J fix)(p{x)dx. 

We write c, Ci, C2, . . . for different constants the meaning of which is 
clear from the context. 

Lemma 2.1. /(see, e.g., [Mu, pp. 57-59])] 

(i) If a; = ayb, where y G 97tn,m; 0. G 6*1/(77,, R), and b G GL(m, R), 
then dx = \a\"^\b\'^dy. 

(ii) If r = qsq', where s G Vm and q G GL{m, R), then dr = \q\'^~^^ds. 

(iii) If r = s~^, s G then r G "Pm and dr = \s\~"^~^ds. 

The Siegel gamma function associated to the cone Vm is defined by 
(2.6) r^(a) = / exp(-tr(r))|r|"-'^rfr, ci = (m + l)/2. 



[Sij . [Gij . |Muj . |FKj . |T] . This integral converges absolutely if and 
only if Re a > d — 1, and can be written as a product of ordinary 
F-functions: 

(2.7) TUa) = vr"(™-i)/^F(«)F(« - i) . . . F(a - ^). 

For n > m, let y„^m, = {f G : v'v = /„} be the Stiefel manifold 

of orthonormal 777-frames in R". We fix the invariant measure dv on 
Vn,m [Mu, p. 70] normalized by 

(2.8) an,m = dv 



v„. 



Tmin/2) 



and denote d^v = a^\^dv. A polar decomposition on 9Jl„,rn is defined 
according to the following lemma; see, e.g., |Mu| pp. 66, 591], |Maj . 

Lemma 2.2. Let x G S!H„^m, n > m. //rank(x) = 777, then 

X = vr^^^, V G Vn,m, r = x'x e Vm, 

and dx = 2-"'\r\^''-"'~^^/^drdv . 

The following statement is new and suggestive. It contains a ma- 
trix generalization of the relevant formula by Smith and Solmon 
Lemma 2.2] corresponding to the case m = 1. 
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Lemma 2.3. Let 1 < k < n — m. Then 

m. 



(2.9) 



f{x)dx 



(^n—k,: 



f{iz)\ztdz. 



Proof. Let / = J^^^f{x)dx, d = {m + l)/2. By Lemma IT^ 

= j dl J f{lVor'/')\rr/'-'dr, Vvo e V^,^. 

SO{n) Vm 

Choose vq = ^0^0, where 


In-k 

By setting ^ = •yC^o, we obtain 
cr, 



Co 



n—k,m- 



' n,m 



j d^i j f{iu,r'/')\r\-/^-''dr ^ ^(3) 



n,n — k 'Pm 



a. 



n,m 



Vn,n-k 0{n-k) Vrr 



<Jr, 



2™cr. 



n—k,m, 



n,n — k ^n — k.m T^m 



^n—kj 



d^i 



f{iz)\ztdz. 



□ 



2.2. Riesz potentials. We recall basic facts from |0R2j and |Ru6j 

related to Riesz potentials of functions of matrix argument. These 
potentials arise in different aspects of analysis |Gej . |Khj . |Stlj . They 
have a number of specific higher rank features and coincide for m = 1 
with classical integrals of Marcel Riesz |Rulj . |SKMj . |St2j . In the 
following, we assume m > 2. The Riesz potential of order a e C of a 
function / G ${^n,m) is defined as analytic continuation of the integral 



(2.10) 



{IV){x) 



/(a;-i/)|y|°-"ci?/, 
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where \y\m = det(?/'?/)^/^ 

(2.11) 7„,r„(a) = ^ ^^^^ , a ^ n-m + l, n-m + 2 



r^((n-a)/2) 

cf . (jl.ip . This integral converges absolutely if and only if Re a > m — 1 
and extends to all a G C as a meromorphic function whose only poles 
are at the points a = n — m + l,n — m + 2, . . . . The order of these 
poles is the same as in Tm{{n ~ Ci)/2). 

Theorem 2.4. // / and (j) are Schwartz functions on dJln,m., then for 
all complex a ^ n — m + 1, n — m + 2, . . . , 



(2.12) 



(r/,0) = (27r)-"-(|y|-°(.F/)(y), {T<t>){y)). 



the expression on each side being understood in the sense of analytic 
continuation. 

This statement is a consequence of the relevant functional equation 
for the corresponding zeta distributions, see |Ru6j . |FKj . 

If a = fc, A; = 0, 1, 2, . . . , m — 1 and /ct^tt, — m + l,^ — m + 2, 
then /"/ is a convolution with a positive measure supported by the 
manifold of all matrices x of rank < k. Combining this fact with the 
case Re a > m — 1, we introduce the Wallach-like set 



(2.13) 
where 



W 



Wi n W2 



1, n — m); 



Wi = {0, 1, 2, ... , fco}, ^0 = min(m 
W2 = {a : Re a>m — 1; ctT^n — m + l,n — m + 2, ...}. 

An analog of ()2.13p is defined in |FKt p. 137] for distributions of 
different type. For a G W„^m, one can write 



(2.14) 



(r/)(a;) = (/*/x„)(a;) 



f{x - y)d^a{y), 



UR2| Theorem 3.14], jRu6| Theorem 5.1], where the measure /Xq, is 
defined by 

1 



(2.15) (/i„,V) 



In- 



[a 



|l/|°X?/)d?/ ifi?ea>m-l. 



Cfc 



dy j ( 7 

0(n) 



y 




(i7 if a = k. 
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Here, k = 1,2, ... ,n 
function, and 



m, ijj is a compactly supported continuous 



(2.16) 



-km 



'km/2 p 
1 r 



n — k 



/r„(i). 



Owing to (|2.12p . for a = 0, fia is the usual delta function, and we set 
I^f = /• Note that the sets of a in both lines of (j2.15j) may overlap. 
In this case we have two different representations of {^a,4')- 

One can use ()2.14|) as a definition of a G W„.m, for arbitrary 
locally integrable function provided the integral f * fia converges abso- 
lutely. 



Theorem 2.5 (jEu6], Section 5.3). Let f e Lp(OT„,„), 1 < p < 

n/ {Re a + m — 1). 

(i) If Re a > m — 1, then 



(2.17) 



exp(— tr(2;'a;)) {I"f){x) dx 



<Ci 



ii) If a = k, k = 1,2, . . . ,n — m, then 

\ii'fKx)\ 



(2.18 



provided 



\Im + X'xl^/^ 



dx < C2 



X > k + max { m — 1 



n + m 



pi 



1 
P 



1 

p' 



1. 



This statement shows that for a G W„ m and / G L^, the definition 
()2.14|1 is meaningful, i.e., {I°'f){x) is finite for almost all x G 9Jl„,,m 
provided 1 < p < n/ {Re a + m — 1). The last equality agrees with the 
classical one 1 < p < n/Rea for m = 1 |St2j which is sharp. We do 
not know whether the restriction p < n/{Rea + m — 1) is necessary if 
m > 1. 

2.3. Radon transforms on the space of matrices. The main ref- 
erences for this subsection are |ORlj . |0R2j . |Pej . We fix positive 
integers k, n, and m, < k < n, and let Vn^n-k be the Stiefel manifold 
of orthonormal {n — A;)-frames in M". For ^ G Vn^n-k and t G DJln-k,m, 
the linear manifold 

(2.19) r = r(e, t) = {x e 97l„,^ : ^'x = t} 

will be called a matrix k-plane in dJln,m- We denote by 1 the set of all 
such planes. Each r G T is an ordinary fcm-dimensional plane in M"™', 
but the set T has measure zero in the manifold of all such planes. 
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Note that r(^,t) = T{^e',et) for all 6 e 0{n - k). We identify 
functions ip{r) on 1 with functions ip{S,, t) on Vn^n-k x ^n-k,m satisfying 
(p{^9', 6t) = ip{^, t) for all 9 G 0{n — k), and supply 1 with the measure 
dr so that 



(2.20) 



fir) dr 



The matrix k-plane Radon transform f{x) —>■ /(r) assigns to a func- 
tion f{x) on Tln^m a collection of integrals of / over all matrix planes 
r G T. Namely, 



fir) 



Precise meaning of this integral is the following: 



(2.21) 



f{9, 



t 



duj. 



where G SO{n) is a rotation satisfying 





(2.22) 



'■n—k 



n,n—k- 



The corresponding dual Radon transform fir) ^{x) assigns to a 
function ^p{t) on T its mean value over all matrix planes r through x: 



V9(r) 



This means that 
(2.23) 



X edJin' 



ip{C,C'x)d4. 



The corresponding duality relation reads 



(2.24) 



f{x)(p{x)dx 



d*C 



Theorem 2.6 (lOES]). 

(i) The Radon transform /(^,t), / G L^{dJtn^m), is finite for almost all 
(^, t) G Vn,n-k X OT„_fc,m ^/ and only if 

(2.25) l<p<p„ = ^±^. 

k + m — 1 
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(ii) Ifip{^,t) is a locally integrable function on the set Vn,n-k x ^n-k,m, 
\ < k < n — m, then the dual Radon transform (p{x) is finite for almost 

The following statement is a matrix generalization of the so-called 
projection-slice theorem. It links together the Fourier transform ()2.4j) 
and the Radon transform (j2.2H) . In the case m = 1, this theorem can 
be found in [Na, p. 11] (for k = n — 1) and [Ke, p. 283] (for any 

< A; < n). 

For y=[yi...ym]^ ^n,m, let L{y) = span(?/i, ...,?/„) be the span of 
the n-vectors yi, ■ ■ ■ ,ym- Suppose that rank(?/) = i. Then dimL(?/) = 

1 < m. 

Theorem 2.7 ([SEIl, ISEl, fnR2] ). Let f E L\mn,m), l<k< n-m. 
If y E 97tn,m; ond C, is an (n — k)- dimensional plane in containing 
L{y), then for any orthonormal frame ^ G Vn,n-k spanning (, there 
exists b G Mn-k,m so that y = ^h. In this case 

(2.26) (^/)(e6) = .F[/(e, ■)](&), hemn-k,n.. 

where T stands for the Fourier transform on Mn-k,m in the {■)-variable. 

Corollary 2.8 ( |UR2 j). The Radon transform / — > / is injective on 
the Schwartz space §(97l„ „i) if and only ii 1 < k < n — m. 

3. Continuous wavelet transforms 

3.1. Some heuristics. Following the philosophy which was described 
in Introduction for the rank-one case, we will introduce continuous 
wavelet transforms on OTn.m associated to the Riesz potential ()2.10|) . 
The heuristic argument presented below shows that these "higher rank" 
wavelet transforms are essentially multiscaled, with the scaling param- 
eter represented by a positive definite matrix, rather then a positive 
number as in the rank-one case. 

We recall the notation ()2.2|) for the invariant measure d^r on Vm and 
start with the following simple observation. 



Lemma 3.1. Let Wq be a symmetric function on Vm satisfying 

(3.1) / I ^!!^°^^|L d^^r < oo and = [ , !^°^^|/^ d^r ^ 0. 

Then for s G Vm, 

(3.2) kl'"-""^ ^ / '"'^^^P . 
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Proof. Using the symmetry (j2.3|) and changing variable a = p^^, d^a = 
d^p, we rewrite (|3.2|) as 

Vm 

It remains to set s^^'^ps^^'^ = r. □ 

According to (j3.2j) . for Rea > m — 1, the Riesz potential (j2.1(J|) is 
represented as 

= [ |a|("-")/^t/,a / f{x-y)wo{a-'/'y'ya-'/')dy 

7n,m(a) J J 

(3.3) = / \<'^d*<^ f f{x-ya''^)w,{y'y)dy. 

7n,m(a) J J 

Vm dJln^m 

The inner integrals in these expressions resemble the wavelet transform 
()1.4p and inspire the following. 

Definition 3.2. Let w{y) = wo{y'y), y G Mn,m, be a radial function 
satisfying certain cancellation conditions (which depend on the con- 
text). Let Waiy) = |a|~"/^w(?/a-i/2)^ ^ ^ ^g^n 

(3.4) (W«/)(x) = J f{x-ya'/^)w{y)dy 

TO 

the continuous wavelet transform of f generated by the wavelet func- 
tion w and the "Pm-valued scaling parameter a. 

Note that the symmetry condition for wq indicated in Lemma 13.11 
plays an auxiliary role. It was imposed only for technical reasons and 
not included in Definition 13.21 In the following, this condition will 
appear on the Fourier transform side for J^w. 

If the function wq in Definition 13. 21 is symmetric, then one can write 
(HOD as 

(3.5) (/"/)(a;) = CnAo^.w) j (W,/)(x)|ar/2 4a, 

Vm 

Rea>m — 1] a n ~ m + l,n — m + 2, . . . , 

Cn,m(tt,'W^) = const. This formula is expected to be true for other 
values of a (e.g., for a = 0) if w obeys certain cancellation conditions. 
Of course, if Rea < m — 1, then the integral on the right-hand side 
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of diverges in general, and must be interpreted in a suitable way 
depending on a class of functions / and a choice of the wavelet w. 

One can replace the function w in ()3.4|) by a finite radial Borel mea- 
sure v on 97l„_m so that 

j ipi'jx) dv{x) = J ip{x) dv{x) 

for all 7 G SO{n) and ijj G Cc(9Jl„^m). If z/ obeys some cancellation (for 
instance, {Ti'){y) = on matrices y of rank < m) we call 

(3.6) {^u,af){x)^{f*va){x)= j fix-ya'/')duiy) 

the wavelet transform of / generated by the wavelet measure u. 
3.2. Calderon's reproducing formula. Denote formally 

(3.7) H^J)i^) = j{'y^u,afKx)d,a, 

where, as above, rf^a = \a\~'^da, d = [m + l)/2. The following state- 
ment justifies fl3.5|) for a = and generalizes the classical Calderon 
reproducing formula (cf. Theorem 1 in |R,u3j ) to functions of matrix 
argument. 

Theorem 3.3. Let u be a radial finite Borel measure on Tln^m so that 
for all y G OKn,m of rank m, we have (jFz/)(?/) = Uo{y'y), where Uo{r) is 
a symmetric function on Vm- Suppose that the integral 

(3.8) c„=mni!l / (^^IM,, 

A<y'y<B} 

(A, B G Vm) is finite. Then for f G L'^{mn,m), 

(3.9) cj{x) = I{iyJ){x) = lim [ (W,,„/)(x)4a. 
Proof. For < e < p < oo, let 

(3.10) 4,(i/,/)(x)= j {W,,af){x)d,a 

elm 
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and assume first that / G flL^. Then, by the generahzed Minkowski 
inequality, Ie,p{^, f) ^ ^ L"^ i and we have 

(3.11) ^[4,p(z/,/)](2/)= j J'[W,,af]{y)d.a. 

elm 

By taking into account that J^[V^u,af]iy) = f) iu) i^^a) {y) , and 
{J^Va){y) = J exp{ti{ia^^'^y'x))du{x) = {J^u){ya^^^), 

we obtain 

(3.12) ^[/,>,/)](y) = My)(^/)(y), 
where 

plm pim 

ke,p{y) = j {^v){ya^^'^)d^a' = j Uo{a^^'^y'ya^^'^) d^a. 

Suppose that rank(?/) = m (the set of all such y has a full measure 
in aH„^m). Since Uq is symmetric, then UQ{a^^'^ra^/'^) = Uo{r^^'^ar^^'^), 
r = y'y, and the change of variable s = r^l'^ar^l'^ ^ yields 

plm pr 



sl„ 



2- f (-Fz/)(z) 

■dz, r = yy. 



n 



n.m'' srKz' zKpr} 

Since ke^piy) is bounded uniformly in e, p, and y, then by the Lebesgue 
theorem on dominated convergence, 

(3.13) ||/,,,(z., /) - cJh = (277)-"™/! (fc,,, - c.)^fh ^ 

as £ — 0, p ^ oo. This proves the statement for f E (1 L"^. A 
standard procedure allows us to extend the result to all f E L"^. We 
recall this argument for convenience of the reader. For any / G L^, we 
have 

plm 

||/.>,/)||2< j ||W,,,/||2rf*a<Q,,||/||2||z^|| =C,,,,,||/||2 

elm 

where \\i>\\ stands for the total variation of c^^p = const, c^^p^u = 
Ce,p 1 1^1 1- Given a small 5 > 0, we choose g E L^f\L'^ so that ||/ — (7112 < 
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6. Since ks,p is uniformly bounded, then ()H.12|) (with / replaced by g) 
implies the uniform estimate 

||/e,p('^,^7)||2 <C||(7||2 <C||(7-/||2 + C||/||2, 

and, therefore, 

\\IeA^J)\\2 < \\IeA''J-g)h+\\IeA^,9)h 

< 5ce,p,u + c\\f-gh + c\\fh 

< 5(q,,,, + c) + c||/||2. 
Assuming 5 — > 0, we obtain ||/e,p(z^, /) II2 < c||/||2. This gives 

\\Ie,pi^J)-'^'^fh < WhA^J - 9)\\2+ \\Ie,pi^^9) - Cj^gh 

+ c^Wg-fh 

< (c + c^)5+ ||/e,p(z/,5() - c^5(||2. 
Owing to (with / replaced by g), 

\\Ie,p{^, f) - Cufh ^ {c + Cu)S as e ^ 0, p -> 00. 
Since 6 is arbitrarily small, we are done. □ 

4. Inversion of Riesz potentials 

We recall that the wavelet transform of a function / on Mn,m is defined 
by 

(W„/)(x)= j f{x-ya'/')w{y)dy={f*wa){x) 

m 

where 

(4.1) Wa{x) = |a|-''/^w(xa"^/2), x G mn,m, a G V^. 

Owing to ()2.12|1 . it is natural to expect, that the inverse of the Riesz 
potential ()2.14|) can be obtained if we formally replace a by —a in 
(ESI); cf. (HISI). This gives 

(4.2) fix) = c„,„(-a, w) I ^^ff ^""^ 4a. 
Below we give this formula precise meaning. 

Theorem 4.1. Let a G W„^m and f & L"^ (1 for some p satisfying 

1 <p < 



Re a + m — 1 

Suppose that w G §(Mn,m) is a radial function such that 

(a) {J^w){y) = uo{y'y), where uo{r) is a symmetric function on V„ 

vanishing identically in a neighborhood of the boundary dVm,; 
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(b) The integral 

,4,3) 1^-^^. 

'-'n.m J \ \in 



B~*oc (J 



\n+OL 



n,m J l^lm 

z'z<B} 

{B G Vm) is finite. Then 

plm 

(4.4) d^{a)f = / d,a = lim / d,a. 

Proof. We observe that the integrand in ()4.3|) has no singularity at 
1^1 = thanks to the assumption (a) above. Moreover, for Re a > 
m — 1, this integral is finite automatically, because 



\y 



dy<c I + y'y\-^''+°'^^^dy < oo, c = c{w 



D3t 5T{ 

see formula (A. 6) in |UR2j . 

To prove the theorem, we set 

pir, 



(4.5) (T» (x) = I ^^^^p^ 4a, < e < p < oo 



el„ 



(if Re a > m — 1 one can assume p = oo), and show that 

(4.6) T,y"/ = ^-i[^,^^^/], 



(4.7) ^,-,^(y) = / 



|n+a 



e(y'y)<z'z<p{y'y)} 

As we have shown this, by the Lebesgue dominated convergence the- 
orem, owing to (j4.3p and the uniform boundedness of ipfp, we obtain 
the desired result: 

||r,y"/ - dUa)fh = (2vr)-™"/l(^°^ - dUa))J'f\\2 - 

as £ ^ 0, p ^ oo. 

We observe that for any f E and w E L} , 

(4.8) T,y"/ = rg, g = T^J. 
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The validity of interchange of integrals follows by Theorem 12 .51 ac- 
cording to which, the integral 1^°"^ [|/| * is finite for almost all 
X because |/| * \wa\ G U'. 

We first prove ()4.6p for / belonging to the Schwartz space § = 
S(9Jl„,,„). By dUD, 

where 

{J-'Wa){y) = |a|~"^^ J exi)(ti{iy'x))w{xa~^^'^)dx 

exp{tT{ia^^'^y' z))w{z)dz = {J^w){ya^^'^). 
TO 

Hence, 

1/2 

iJ^g)iy) = KM{Tf ){y), KM = j ^^^^^^ d^a. 

elm 

By taking into account that {Tw){y) = Uo{y'y) where Uo{r) is symmet- 
ric, we obtain 



he,p{y) = I rrp^ 







I 




elm. 




2m 


\y'y 



I 



{zeM e{y'y)<z'z<p{y'y)} 



dz 



(4.9) = lyC^fM 

(see the argument in the proof of Theorem IH.Hp . This gives 

(4.10) iJ^g){y) = \y\l^:MiJ^f)iy). 

Since w G § and Uq is supported away from the boundary dVm, it 
follows that {J-'g){y) = he^p{y){J-'w){y) G §, and therefore, g & §. 
Hence, by (14. Sp . ()2.12j) . and ()4.1()j) . for any compactly supported C°° 
function 0, we have 

{T,y^f, 0) = {rg, 0) 

= (27r)-"™(^r^, JP-(/.) 

= (27r) — (|y|-"(^^7)(y), (^0)(y)) 

= (27r)— (^,",^(?/)(^/)(|/), (^0)(y)). 
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Thus, by the Parseval equahty, 

(4.11) (T,y"/, 0) = 0) 

(note that V'"p(?/)(-^/)(l/) ^ L"^ in view of the boundedness of 'ip^^piy))- 
Since J-'~^[ip^ ^J-'f] G and T^^I'^f = I°'g is a locally integrable func- 
tion (see ()2.17|) and ()2.18|) ). then ()4.1ip implies the pointwise equality 
for any / G S. 

To complete the proof, it remains to extend (j4.6|) to all / G fl L^. 
Following Theorem 12.51 we introduce the weighted space 

X = {if : WfWx = j \(f{x)\ ijj{x) dx < oo} 
TO 

where to'(x) = exp(— tr(x'x)) iiRea > m — 1, and a;(x) = \Im + x'x\~'^^'^ 
if a = k, k = 1,2, . . . ,n — m; see (|2.18p . It may happened that these 
domains of a overlap, but this is not important. By Holder's inequality, 
I Iv^l U < c I Iv?! I2. Since bv Theorem 12.51 

WT^J^Wx = \\I"g\\x <c\\g\\, = c\\T^J\\, < c.JJX, 

and 

operators 

are bounded. This remark allows us to extend ()4.6|) to all f & L'^ H 
W by taking into account that there is a sequence {fj} C § such 
that the quantities ||/ — fj\\p and ||/ — fj\\2 tend to as j — 00 
simultaneously. Such a sequence can be explicitly constructed using 
the standard "averaging-truncating" procedure. □ 

5. Continuous ridgelet transforms and inversion of the 

Radon transform 

5.1. Intertwining operators. Given a sufficiently good function w 
on Mn-km^ consider the intertwining operator 



(5.1) {Wm,t) = j f{x)w{t-i'x)dx 

which transforms a function / on Ttn,m into a function Wf on the 
"cylinder" Vn^n-k x ^n-k,m- The corresponding dual operator is defined 
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by 

(5.2) (W^V)(a;) = J j ip{i,t)w{i'x-t)dt, 

SO that 



(5.3) J fix)iW*if)ix)dx= J J ipi^,t){Wm,t)dt 

(at least, formally). We shall see that Pm-scaled versions of W and W* 
can be regarded as matrix modifications of continuous /c-plane ridgelet 
transforms (see |('aj . |R,u5j . and references therein), and used for ex- 
plicit and approximate inversion of the Radon transform 1)2.211) . We 
start with some preparations. 

Lemma 5.1. Given a function f{$,,t) on Vn,n-k x ^n-k,m, l^t 



(5.4) (Wom,t)= J ip{^,z)w{t-z)dz 
be a convolution in the t-variahle. Then 

(5.5) {wm,t) = {Wom,t), 



(5.6) 



{w*v){x) = {w,^y{x) 



provided that either side of the corresponding equality is finite for f , ip, 
and w replaced by \ f\, \ip\, and \w\, respectively. 

Proof. The equality ()5.6|) follows immediately from ()5.2|) . To prove 
(15. 5p . we choose a rotation e SO{n) satisfying 







ri'n—k- 



The change of variable x = g^y in ()5.1|) gives 

{Wf){^,t)= [ f{g^y)w{t-i',y)dy. 



By setting 



UJ 

z 
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SO that = z, owing to (j2.2H) . we obtain 



w{t - z) dz J f [g^ 



UJ 

z 



f{^,z)w{t-z) dz=iWom,t). 



□ 



Lemma 5.2. Let f{x) and w{z) be integrable functions on ^n,m one? 
SWn-fc,m; respectively. Then {W* f){x) is a locally integrable function on 
^n,m which belongs to S'{Mn.m) (irid satisfies 



(5.7) J iW*f)ix)(l){x)dx= J fix)iW*(l))ix)dx, 0eS(OT„,, 
Proof. We have 



l.h.s 



d^ 



(5.i 



f 2nv 



f{u){wm,t)dt 

f{U){Wom,t)dt 
fix) {Wo(j)y{x)dx 



r.h.s. 



These calculations are well justified and all statements of the lemma 
become clear, owing to the following estimate of the expression (j5.8|) : 



(5.9) 



< 



4e J \fi^,z)\dz J \(l)i^,t)wit-z)\dt 
oolklli j d^ j \f{i,z)\dz 



l«^llill/lli- 



□ 
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In the sequel, it is convenient to use different notations for the Fourier 
transform on Mn^m and Mn-k,m- For the first one we write JF as before, 
and the second will be denoted by JF. 

Lemma 5.3. If w G L^{Tln-k,m) and G S(9Jl„^m); then 
(5.10) {W*^){x) = (27r)('=-")"^ j 4e 

X j ey.]i{—ii{ix'^z)){Tw){z){T(j)){iz)dz. 

Proof. Since w G L^($H„_fc,m) and G S(9Jl„,m), the convolution 



{W,m.t)= j <P{i.t-z)w{z)dz 

has the Fourier transform (in the t- variable) belonging to L}{Mn-k,m)- 
Hence (see, e.g., |SW| p. 11]) one can write 



{Wom.t) = i2nf-''^"' j expi-tTitt'z))[:FiWom,-mdz 

= (27r)('=-")™ j exp{-tT{it'z)){:Fw){z)[:F^i^,-)]iz)dz. 

By the projection-slice theorem (see ()2.26p ). 

{Wo m,t) = (27r)(^-")™ I exp{-tT{zt'z)){:Fw){z){J'mz)dz. 

^jVfi — k,,m 

This proves the statement. □ 

5.2. Continuous ridgelet transforms. Let w{z) be a sufficiently 
good function on dRn-k,m^ \ < k < n — m. We consider the Vm- 
scaled version of w defined by Wa{z) = \a\^^~'^''^'^w{za~^/'^), a G Vm, 
and introduce the following dual pair of intertwining operators 

(5.11) iWafmt) = j fix) Wait- ex) dx, 



(5.12) (iy»(x) = 

V„,r^-k Tl 



j 4^ j ^ii,t)Waii'x-t)dt. 
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If w{z) oscillates in a certain sense we call (j5.11|) the continuous 
ridgelet transform of /, and (j5.12j) the dual continuous ridgelet trans- 
form of !f. 

Operators ()5.1H) and ()5.12|) generalize usual fc-plane ridgelet trans- 
forms |Caj . [Ru5j to the higher rank case m > 1. The function x 
w{^'x — t) is constant on each matrix plane r = {x G DJln,m '■ C'x = t} 
and represents a "plane wave" . 

The following definition will be useful in the sequel. 

Definition 5.4. A function w{z) on Mn-k,m is called an admissible 
wavelet function if it obeys the following conditions: 

(i) w{z) is radial, i.e., w{z) = wo{z'z), and belongs to L^(Tln-k,m)- 

(ii) The Fourier transform of w has the form {J-'w){y) = Uo{y'y), where 
Uq satisfies the symmetry condition ()2.3|) . 

(iii) The integral 

(5,13) c„ ^ Hli^ / 

^n,m J Isim 

= hm / ^—^dC 

{f6ajl„_fe,™:A<C'C<B} 

{A, B e Pm) is finite. 

5.3. Inversion of the Radon transform. 



5.3.1. Discussion of the problem. There exist different approaches to 
inversion of the Radon transform (|2.21|) : see |UR2j . The consideration 
below sheds new light on this problem and provides essential progress. 
To explain our strategy, we use intertwining fractional integrals P"/ 

and P°'(p of the Semyanistyi type, which link together the Radon trans- 
form f{x) f{^,t), the dual Radon transform ip{^,t) 'fix), and 
Riesz potentials. Namely, we define 

(5.14) PV = rf, PV=(/»\ 

aGC, OT^ra — A; — m-|-l, ra — fc — m-|-2,.... 

Here, 1 < k < n — m and /" denotes the Riesz potential on Mn-k,m. 
acting in the t- variable. Operators ()5.14p were introduced in |0R2j . If 
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Rea > m — 1, they are represented as absolutely convergent integrals 

(5.15) (p°/)(e,t) = — / /(x)ie'x-t|°+^-"rfx, 

(5.16) (P»(x) = [ d4 L{i,t)\i'x-tC'--dt, 

where 7n-fc,m(«) is the normalizing constant in the definition of the 
Riesz potential on Mn-k,m', cf. flLlip . ()2.11|) . 

The following statement determines our way of thinking. 

Theorem 5.5 ([OaS], Section 5.3). Let 1 < k < n - m, a e Wn,m- 
Suppose that 

Ti 

Rea + k + m — 1 

Then 

(5.17) (p"/)(a;)=c„,fc,„(/"+V)(a;), 



(5.18) c„,,,^ = 2^-vrW2r„ (^^j /r„ . 
In particular, for a = 0, 

(5.19) ifyix)=Cn,kAl'f)i^) 
(the generalized Fuglede formula). 

Formula ()5.17p paves two ways to the inversion of the Radon trans- 
form. Following the first one, we use ()5.17|) as it is and invert the Riesz 
potential 1°^+^ f by choosing a as we wish. For instance, one can set 
a = and apply (|5.19|) . This program can be realized using results of 
the previous section. The second way is to set (formally) a = —k in 
()5.17|) . This gives 

(5.20) Cn,k,mfix) = ip-'f)ix) = ii-'fnx), 

and we have to find "good" representation for the inverse of the Riesz 
potential applied to /(^, t) in the t-variable. In the first case, we just 
apply the left inverse operator to J°+'^. In the second one, we do not 
know in advance whether /(^, •) lies in the range of the Riesz potential 
I^. To circumvent this difficulty, we make use of continuous ridgelet 
transforms. 

Below we consider both approaches. 
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5.3.2. The first method. We utilize the generalized Fuglede formula 
(j5.19|) and invert the Riesz potential according to Theorem 14.11 This 
gives the following result for the Radon transform. 

Theorem 5.6. Let 1 < k < n — m, 

(5.21) feL'nL^, 1<P<7— 7- 

k + m — 1 

Suppose that is the continuous wavelet transform ^3.4\ ) generated 
by the wavelet w satisfying conditions of Theorem \4.1\ Then the Radon 
transform f{x) f{^,t) can be inverted by the formula 

f 'y^af f Waf 

(5.22) d^{k) f = Cn,k,m J j^j^ d^a = Cn,k,m J |^p^c?*« 



where 

5.3.3. The second method. By ()4.4|) and ()5.20|) . it is natural to expect, 
that the Radon transform can be inverted as 

... f,, f {f{^,-)*^am'x) , f {w:f){x) ^ 

f{x)= jd^ij d,a= J ^^^d,a 

(up to a constant multiple) where W* is the dual ridgelet transform 
(|5.12p . Below we justify this formula. 

Theorem 5.7. Let f G L^{mn,ni) n L'^{mn,m), I < k < n - m. If w 

is an admissible wavelet function (see Definition \5.4\ l, then the Radon 
transform f{x) f{C,,t) can be inverted by the formula 

f W*f ) f W*f 

(5.23) c^f = J d,a = lim J j-^ d,a, 

where defined by \5.1'Jji . and 

(5.24) {W:f){x)= [ rf.e / fi^,t)wa{ex-t)dt. 



Proof. Consider the truncated integral 

r w*f 

(5.25) IeJ= J j^d.o., 0<£<p<oo. 

elm 
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Owing to for any test function G §, the expression {Ie,pf, 0) is 

finite when /, w, and are replaced by |/|, \w\ and |0|, respectively. 
Hence, we can change the order of integration, and ()5.7|) yields 



(5.26) {I, J, 0) = (/, 



pi 



where le^p has the same meaning as in (|5.25|) but with w replaced by 
its complex conjugate w. Let us show that 

(5.27) = 

where m^^p{y) = me,p{y'y), 

,5.28, ^ / (Mo^c. 



ICIS. 

{C69Jl„_fc,„i:er-<C'C<P'-} 

r G Vm- Suppose that the Fourier transform of w has the form 
{^w){Q = Mo (CO where Uq obeys the symmetry condition (|2.3|) . By 
(with w replaced by w^), we have 

(4»(x) = (27r)('=-")- J d4 y exp(-tr(^x'ez))A:,",(^)(^0)(^^)dz, 
where 

Owing to the symmetry ()2.3p . we have Mo(a^^^ra^/^) = Mo(r^/^ar^/^) 
(without loss of generality, one can assume rank(z) = m). Then by 
changing variable r^/^ar^/^ = s and making use of Lemma 12. 2t we 
obtain 

pim pr 



sir, 



^n—k,m J \>>\m 

{CG97I„_fc,„:£r<C'C<pr} 



Replacing ^ by — ^ and using the equality {fiw){—Q = (^w){(), we 
have 
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Hence 



(/e,p0)(x)=c j d4 j exp{-tT{ix'Cz))\z\^ms,p{z'z){J^(f)){Cz)dz, 
where c = (27r)~"'™'(j„^m/c"n-fc,m- By Lemma ESI 



= (27r)-"'" j exp{-tT{zx'y))m,^,{y'y) {:Fct>){y) dy, 

and (jKTTfjl follows. 

The rest of the proof is standard. Since f & and m^^p is uniformly 
bounded, then by the Parseval equality, 

(5.29) {hj, 0) = (^-^m,,,^/, 0), V0 G S(M„,^). 

By Lemma 123 Ie,pf is a locally integrable function. Since T'^m^^pT ] 
is locally integrable too, then ()5.29|) implies a pointwise equality /e,p/ = 
T'^rrie^pT f ^ and we have 

= (27r)-"-/2||(m,,,-c^)^/||2^0 
as £ — > 0, p ^ oo. □ 

5.4. Reproducing formula for the ridgelet transform. Given two 
functions u{z) and v{z) on "^n-k.m-, we set 

and consider the corresponding ridgelet transforms 

(5.30) (f/a/)(e,t) = j f{x)Ua{t-Cx)dx, 

(5.31) {V:^){x) = j d4 j ^{^,t)va{^'x-z)dz, 

cf. dnH), (EI21). 

Theorem 5.8. Lei u and v be integrable radial functions on Mn-k,m, 
1 < k < n — m, such that their convolution u * v is admissible (see 
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Definition \5.4]) - Let 



B^!L ^n,m J ICIm 

{Ce9Jl„_fc,„:A<^'^<B} 

fA, B e Vm). Then for f e L\mn,m) n L^{mn,m), 

f V*U f ) f V*U f 

(5.33) Cu,J = J d,a = lim j d,a. 

I tn ^ fcJm 

Proof. Let us show that V*Uaf coincides with the dual ridgelet trans- 
form W*f (see (|5.24p ) generated by the function w = u * v. We have 

iW:f){x) = J jf{^,z)dzj Ua{t)Va{i'x-Z-t)dt 

= j d,i jf{i,z)dzj Ua{Q-z)Va{i'x-C)dC 

d*i j Va{i'x - C)dC j UaiC - z)f{i, z) dz 

,n — k 9Jt — k,ni ^jytji — k,ni 

j 4e jvaii'x-Qdcj fix) UaiC -^'X)dx 

= iv:uj)ix). 

Now the result follows by Theorem 15 .71 □ 
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